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In this paper we describe the assembly and restric- 
tion map of a 1.05-iVtb cosmid contig spanning the can- 
didate region for familial Mediterranean fever (FMF), 
a recessively inherited disorder of inflammation local- 
ized to I6pl3.3. Using a combination of cosmid walking 
and screening for PI, PAC, BAG. and YAC clones, we 
have generated a contig of genomic clones spanning 
•1050 kb that contains the FMF critical region. The 
map consists of 179 cosmid, 15 Pi, 10 PAC, M BAG, and 
17 YAG clones, anchored by 27 STS markers. Eight ad- 
ditional STSs have been deveioped from the -700 kb 
immediately centromeric to this genomic region. Five 
of the 35 STSs are microsatellites that have not been 
previously reported. Notl and EcoRI mapping of the 
overlapping cosmids, hybridization of restriction frag- 
ments from cosmids to one another, and STS analyses 
have been used to validate the assembly of the contig. 
Our contig totally subsumes the 250-kb interval re- 
cently reported, by founder haplotype analysis, to con- 
tain the FMF gene- Thus, our high-resolution clone 
map provides an ideal resource for transcriptional 
mapping toward the eventual identification of this dis- 
ease gene. ^m^ Academic Press 
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INTRODUCTION 

Familial Mediterranean fever (FMF) is a recessively 
inherited infiamnnatory disorder occurring primarily in 
non-Ashkenazi Jews, Armenians, Turks, and Arabs 
(reviewed in Kastner, 1996). The FMF gene frequency 
among these populations is very high, with estimates 
based on family studies as high as l;iO among some 
non-A.shkenazi Jewish populations (Daniels et a/.. 
1995) and 1 in 14 among Armenians in Los Angeles 
(Rogersei a/., 1989t. The mam clinical features include 
intermittent attacks of fever with abdominal pain, 
pleurisy, and arthritis of one or more joints; patients 
are normal between attacks. Patients may also develop 
systemic amyloidosis, leading to renal failure and 
death. In the absence of a known biochemical abnor- 
mality associated with this disease, our laboratories 
have established a consortium to identify the gene by 
positional cloning. 

Following a genome-wide linkage analysis in a panel 
of non-Ashkenazi Jewish families, we localized the 
gene for FMF (designated MEFV) to the short arm of 
chromosome 16 (Pras et al., 1992). The FMF suscepti- 
bility gene maps to the same chromosomal region in 
Armenians (Shohat et aL, 1992), Arabs (Pras et al., 
1994), and Turks (Ozen et aL, 1996). Subsequent ge- 
netic analyses placed MEFV between D16S94 and 
D16S80, a genetic interval of —9 cM (Aksentijevich et 
aL, 1993a; Fischel-Ghodsian a/., 1993). By the outset 
of the present physical mapping studies, recombinants 



83 



083a-7.=>43/97 $25 DO 
CopyriRht 1997 by .Academic Pre«;s 
All rights of reproduction in any form rp.served. 



84 SOOD ET AL. 

TABLE 1 

STSs from the FMF Candidate Region'* 

Size An: 





STS name 


Location 


Primer sequences^ i'5'~3'> 


(bpi 


te:: 




s;iV0G9 I 


C310G9 


GTGCTCTCCCAACCTCTCAG 


121 










AGTGAGCACCTGGGCTTG 






Dif;S3:i66 


s302H7_l 


C302H7 


CGCAATTAACGCTCACTAAAGG 


1 1 1 










GGGACTC'ITCCAGGAGATCC 






DiryS2S-}4 


?54G6_ i 


CS4G6 


TCT0GTGGCCCTACGAA1TGAG 


388 










G G G TG G G ACriTAG A ATCC G G T 








s5'lG6 2 


C54G6 


GTCGTTATGCCCGCTGAGTATG 


123 










CCTAACCTCTCTGGGCCTGTTT 










o54Gr> 


AAGTCCCAAAAGTCACACCG 


2fi3 










GGAGTTATCAGTGCAGATTACGG 








C28iCA)/C28£GT.' 


c 3 63 09 


CCCTCCTTGAATTACTTGAACACGGGA 


194-216 




'.niGS.wiOf 




0 8 kb 


TG G G CrC C C A G G C AG CTCTG A A AO TO CnTT 






iAC)„ 
D16S.WS2 f.AC)„ 


S49U-1 


<:49B4 9 kb 


TGCGC/VAATAACGGTGACACTC 


118- 156 










TCTCCTCTTATCCTGTACCTCTGCC 








s.'134Bl2_l 


C334BI2 T7 


G TG AG CC ATC ATGCC C AG AC; G 


101 








end 


GTTGCACTGAGCCAAGACCAC 






DI^S3369 


S406E2 1 


C406E2 T3 


A'rrCTGCTCATATGTCACCC 


IIS 








end 


GTCGTGGO/\AACACGTATAT 






Dt6S3370 (AC)., 


S57K1 1 


colt A T7 


CTGGCTAGAGAACCACGTAAAC 


141 -t53 








end 


ACTTTGCTAGTGAGGTGATTGG 






Dl6S26t7' 


CHLC.ATA4 IE08 


C414G4 3,5 


TGTCCACTAGGTAAAGGCTO 


93-117 




(ATA), 




kh 


TGGGCAACAAGAATGAAATT 






DI6S3371 


S4UG4 1 


c4t4G4 T7 


C'TAGATAAGGCCACACTGACTGTGr 


210 










CKiCCATGAATGAGTATCTTAATAT 






D16S3372 


s442B'2 I 


C442R2 T3 


ATTCAACACCCTTCTTCTGGC:CCA 


181 










GAGCTAGTGAGCTGCATGAGCTGT 






DI0S3373 lACu 


.s374H9_l 


C374H9 4.8 


TrnTAAAAGCAGACTCTGCCC 










kb 


TCTGCATAAATGTTATGCCTGC 






D 16^^3275 (AC),. 


AFMer;i4 


C360H6 13.3 


CAAAGCCCTAAAGTAGCAGT 


198- Ki4 








kb 


GGGTTTGGAGATTTTTGTAA 






Ot^S:)374 


s27:iL24 1 


('AC 273 L24 


AACGG/VACACAGGAAACAfTGO 


210 








SP6 end 


A G AG C AAAC CC: AATG C GGTG 






ni6S337r, 


s2:nf;7 I 


iy23lE7 


C A C CGTG Cr rr AAC TTG AAT 


142 








Right end 


TGAAAAATACiCAnTGAGf\ATAA 






n!/iS33 7/U..\C>. 




c38ril)9 4. 1 


T'nAATCCTGGGACATTCTTGG 


197 -227 








kb 


ATC AGGAA( ;(U;GT( If ;G AGr 






[)!h-S3377 


s3(;()Al I 


c3BnAl T3 


gtactgc(.;gtcattg(:tca(.; 


if^3 








end 


GTTTGCTGTTCITGGCCG 






l.)!fiS3378 


sO()4fiC8 1 


VAC O046C8 


(^ACATTTTCAGOfJATTGAGAATTA 


156 








SP6 end 


TGTGAGCCACT(XTCX:TGGGCATC 






D}fiS3379 


S441H9. 1 


c44in9 T7 


GAGACAGGGTCTCGCTATGT 


221 








end 


rTCnTGTCCACTACTTCJGT 







within our panel of 65 families defined the candidate 
interval as te\-'Dl6S246-MEFV~D16S2622, a genetic 
distance of — 1 cM (Sood et aL, 1996). Taking advantage 
of the founder effect observed in the Moroccan Jewish 
population (Aksentijevich et al., 1993b), and applying 
Luria-Delbrlick formulas and simulations based on a 
Poisson branching process, we estimated the gene to 
be within 0.305 cM o(D16S246 (Levy et aL, 1996). 

D16S246 lies at the centromeric boundary of the ge- 
nomic interval from which the polycystic kidney dis- 
ease gene (PKDl) was cloned (Harris et aL, 1994), while 
D16S2622 is derived from a cosmid (cRT70) that is 
near the telomeric boundary of the region in which the 
Rubinstein --Taybi (RSTS) gene was found (Petrij et 
ciL, 1995). Consequently, the interval between these 
markers had not been mapped in the course of these 
two positional cloning projects. Moreover, this region 
also lacked clonal representation in the integrated 
physical map of chromosome 16 developed by Doggett 
et aL (1995). By screening a combination of YAC, BAG, 



PAC, PI, and cosmid libraries, we developed a hi 
redundant, internally consistent physical map of 
~1050-kb region spanning the MEFV candidate it 
val. Microsatellites identified in the process of 
structing this map have permitted further narro\ 
of the interval (J. E. Balow et aL, in preparation). 

A second FMF consortium comprising several Frt 
groups recently published data localizing MEFV tc 
250-kb YAC clone 26Fe7 (iy23lE7; French FMF ( 
sortium, 1996). Our own genetic analysis and the ' 
in the present report are in agreement with this ; 
ing, but provide a much higher resolution map of 
region. By placing our own and the French consortiv 
genetic markers on the EcoRl restriction map, we 1 
created an integrated physical/genetic map of 
MEFV candidate interval, 

MATERIALS AND METHODS 

YAC libraries and ficreening. YAC clones were isolated by 
based screening ofseveral different total human genomic YAC li 
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TAIiLE l—Contmui>xl 

Locus STS name L,Dcatinn [*rimei- suquences'' (f)' ~,T i (bp) temp i^T; 



D16S475 


a58H 


iy58H4 


TCfCAtiCvAAl AH.XrAl(,'OA(_r t 1 ^iVj 




45 


(AGlJ AAAGi, 




[.eft end 


GOAGAGGAAGAGCGATGGGAGTAT 






D 1633380 


93o8BlO_ i 


T7 (Mid 


GGC AC AGAGG 1 1 G [ GG t GG( i 
G C A CTTCGTC CO CCTCTC; G( ; C 






!) lasses ! 


sfn'21 ] ^ 1 


rRT21 1 


gcaacagagcaagactctgc 
gtgaggtgacgcaactagagg 


129 


■5 ? 


D16S3382 1 AC^ 


s400C4_l 


C-100C4 


ACTTCTAAGCTCCTGACATGC'C 
ATGATACTTCTGCTGTAAGCGG 






l)!6S3383 


s.'iOTEfi I 


c307Efi 


tgaaaaagctgctcacatc 


1,15 


""^ 


l)HiS2622 


-sRTTO 1" 


rirrvo 


TGACTGTAGCTTGGGTGAv^GG 


156-168 










cctctcgagaggagaactgg 






ni(iS33S'i 


sRTa_ 1 


cRTS 


AATGAAGTAAAAAGTGGrnTGGG 
CTAGAGGGAAGCACrrTGTGG 


172 








cKT-t 


ACGGAAGCTAriTGGfK^C 
AAAAC GOT A A AA CTG A G AGG AG G 


159 




D!6S290e 


sKTlJ 


cIlTl 


ATTTAAGTGAri'GGCAGCAGGG 
GC/^AGAAACAACAGCCTGGA,^G 


399 




D16S290T 


sRTl 2 


cRTl 


T GTTTGG AG C CTG TGCTGITTG 
CCACTCTGAAAAG^^ACAGCCGG 


ion 




D 16833 8ft 


sRT5,1 1 


cRTf>3 


C G CTG GTTC ATGGTAA C ATG 
AAGrrGAACAGCTGTGTCTAGTGG 


222 




D16S3387 


sRT99_ I 


cRT99 


CTCAGCTCAGGATGCTCTCC 
A GTC CTGT A G GGT AG A GG AG GG 


14t 




DieS338H fAC)„ 


sN2_t^508 


cN2 


ACAACCCTGCTTACACCCTG 
GGGA.V(\TTCCATCTCCACAy\ 


143- 171 


52 


D16S3389 


sRTlfi3_l 


cRTlF,3 


CACAGGCACAAGCACGAG 


Uo 





" Listed in their genomic order, telomet-ic to centromeric. 
In cases where there are multiple primer sets for a given locus, we have listed the primers and conditions used in the present experiments 

' D16S468 is the same locus as 01633070 (AFMa35.3yh U. Both map to the same ^coRI fragment, and the C28(CA)/C28(GTl forw^-ird 
primer is nt 120-145 of GSDB Accession No. Z5r;0L3 (the sequence from which D16S3070 primers were derived). The reverse primer of 
C28(CA)/C28i'GT) is outside of Z53013, 

''Primers shown here were developed by screening c49B4 for niicrosatellites. We subsequently found that s49B4_l and AFMbO70yc5?i 
{Dl6S:i082) map to the same Ecolii fragment. Comparifion of 54984 1 sequence with that of GSDB Accession No. Z5.3260 ithe seqLiencv 
from which AFMbOTOygSa primers were derived" showed a 96% identity over 272 bp that contained the (AC)n repeat. Therefore, both 
of primers identify 

' Equivaletit to ni6S2(U8 tCHLC ATA42A()«.P;i4342l, The Dl6S2(il8 forward primer is nt 2 -21 of GSDB Accession No. Gl(J25o rht 
sequence from which D16S2617 primers were derived:'. The reverse primer o{ DIHS2618 is the reverse complement of nt litfi- 116 of (; l<.»2rS3 
Primers shown in thi.-s table wei-e derived from end sequence of iy58I-I4. Later, we found that GSDB sequence Lifi287 ifrotn wfuch. 
D16S475 (UT581) primei'S were derived) is included within this end sequence. The r)8H4 forward primer shown here is the i-e\er=-j 
complement of nt 382 -404 of L 16287. The 58 H4 reverse primer is outside the limits of L 16287. The allele sizes listed here were obtamed 
by digesting the amplification products with Hindll. 

''5RT70 1 and CHLC.GATA73G0o are both associated with D16S2622. The primers shown here were developed by screening cRTTO for 
microsatellites. On subse{iuent screening of other new marker? reported from tl^e region, we found the CHLC.GATA73G05 primers w.thm 
the sRT70 1 amplimer sequence. 



ies. These included the CEPH megaYAC librar)- (Bellanne-Chantelot 
et a!., 1992; clone names prefixed with AfyK the Washington Univer- 
sity YAC A and B libraries fBrownstein ci ai, 1989; clone names 
prefixed with Wy), and the ICI YAC library (Anand e.t ai, 1990; clone 
names prefixed with iy). Agarose plugs of yeast cells containing total 
yeast and YAC DNA were prepared (Gemmill ct al., 19951 and used 
to estimate the size of YACs by pulsed field gel electrophoresis. DNA 
minipreps of YAC colonies were performed using Puregene kits (Cen- 
tra) and were used to confirm STS content mapping and as templates 
for the generation of inter-A/« PGR products (Liu et al.. 1993) for 
hybridization probes. 

BAC, PAC, and PI libraries and screening. BAC clones (Shizuya 
et al., 1992) were identified by hybridization of cosmid EcoRI frag- 
ments to high-density gridded human genomic BAC filters purchased 
from Genome Systems. PAC clones (loannou et al., 1994) were identi- 
fied by PCR-based screening of a pooled human genomic PAC library 
(Genome Systems) with STSs from the MEFV interval. PI clones 
(Shepherd eC al., 1994) were identified both by PCR-based screening 
(clones prefixed with A) and by hybridization of cosmid restriction 
fragments to high-density gridded human genomic PI filters (clones 
prefixed with B). 



Cosmid libraries and screening. Cosmid C28, which contains the 
microsatellite marker C28(CAi/C28(GT) (Di65^6S). wafj kindly pro- 
vided by Dr. K. Hayashi (Kyushu University, Tokyo). Cosmids cRT4, 
cRT8, cRT70, cRT194, cRTl97, and cRT21 1 from the RSTS candidate 
region (Petrij et ai, 1995) were kindly provided by Drs. Rachel Giles 
and Martijn Breuning (Leiden University, The Netherlands), while 
cosmid cN2 was provided by Dr. Anna Maria Frischauf (Imperial 
Cancer Research Fund, London). All other cosmids were identified 
from Los Alamos National Laboratory flow-sorted chromosome- 16 
cosmid hbraries by hybridization to high-density gridded mem- 
branes. These included ~4OO0 fingerprinted cosmids (Stallings ct aL, 
1990, 1992) arrayed on 3 membranes and —15,000 cosmids ! Long- 
mire el al., 1993) arrayed on 10 membranes (library 16-2). Through- 
out this paper individual cosmid names are prefixed with c. and 
cosmid contigs are prefixed with C. 

Inter-A/« PGR products of YAC, PAC, BAC, and Pi clones, and 
rihoprobes generated from the ends of cosmid clones, were used as 
hybridization probes after blocking for repetitive sequences using 
1000-fold excess of Cotl DNA (Life Technologies). Stepwise cosmid 
walking was performed by generating a T3 or T7 riboprobe from the 
farthest extending end of the cosmid clone. One microgram of HaeUl- 
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digusted cop mid DNA was ui^c-d Lo gotiorntc iihnprf>tir' ( Maxiscript in 
vitro trnn.striplinn kit; Artiliioni in tlio proscnce of 50 //Ci of lf^-^''P|- 
tJTP fspecific acliviLy 800 Ci/nitnol). DNA from all positive clones 
was purified from l.'^O-ml cultures us j tig a coninierciatly avaihiblo kit 
fQiagen). One microgram of cosniid DNA was digested with EcnR\, 
separated on 0.9,9r agarnse gols. and blotted onto Gene Screen Plus 
(Dupont) membranes. The overlap of cof!intds (indicated hy restric- 
tion digest? was verified by hybridization EcoKX fragments of 
the starting cosmid. 

Restriction niapping. Detailed EcoWl and Not I restriction maps 
spanning the MEFV region were detet mincd by end-labeled re- 
striction mappinf? of all cosmids with the olign mapping method 
described in Evans et al. { 19H9), Briefly. 4 /'g of cosmid DNA was 
digeKted with Not I, and an aliciuot ( 1 //gi was saved to detect inter- 
nal Noil sites. 'I'o the roniainder of iVoCl-digcstod DNA, 5 j/1 of a 
1/100 dilution ofEcoRI (10 U'V/l) m 1 <■ buffer was added and incu- 
bated at .^VC', A li (pilots of i //g each were removed after I and 5 
min of incultation, and digestion of the remaining DNA was carried 
to completion by adding 1 ^1 nf^raRI. in addition, 1 ^g of cosmid 
DNA was digested with EroKI alone to Incalixe the Nofl sites 
within £coRI frngments. All samples were separated on 0-8^'' aga- 
r<ise gels, blotted onto Gene Screen Plus memlM-anes, and hybrid- 
ized with ■''■^P- labeled T3 and T7 oligonucleotides sequentiaUy. The 
DNA Analysis Marker System ^ Life Technologies 1 was used to esti- 
mate the Size of the partial fragments hybridizing to the T3 and 
T7 oligonucleotides. 

Probe labeling and hybridizations. Oligonucleotide probes wei-e 
end-labeled using T4 polynucleotide kmase and j > - ''P}ATP and were 
passed through G 25 -Sep hade. x spin columns (5 Prime -3 Prime) to 
remove the unincorporated nucleotides. IrUer-A/:< PGR probes and 
EcoRX fragments from cosmids were labeled using a random-priming 
kit fStratagene) and In-- '"fMdC'TP. All pi-obes eKcept oligonucleotides 
were blocked with an excess of r„/l DNA prior to hybridizations. All 
Southern blot hybridizations wore performed in Rapid- Myb buffer 
(Amersham Life Science!. The blot.^ were prehybridizod for L'5-30 
min, hybridi?,{'d for 1-2 h. atui washed twice at room temperature 
in 2x SSC, 0.1 SDS folltnved by washing at firVG in 0. 1 x SSC. 0. 1% 
SDS. For gridded cogniid libraries, hybridizations were performed 
overnight in 0.5 A/ sodium phf>.';phate * pf I 7.2 i. I'^.y SDS, 1 ruM EDTA, 
followed by wasiies as desctdH.-d above 

Cri'iicratinri nfSTSs. • i ) iMid-clrnie STSs: T.'! or T7 ends of cosmids 
were sequenced to generate 9 number of STSs for screotiing of large- 
insert clone librarie.s. Etid fragments of YAG clones [y'tSHA (7Dli4) 
and iy231L7 (26Fe7l and of PAG 2731/24 were obtained by bubble 
PGR using linkers complementary to IfacUl and flsal sites (Kilcy ci 
al., 1990). Oligonucleotide primers were selected with the aid of the 
computer program PRIMER i Whitehead Institute) or ASTS (Rappa- 
portet al.. 1994). I'C-R conditions for alt STSs are described in Table 
1. Chromosome- Ifi-spocific breakpoint hybrids CY18, CY14, CY18«, 
23HA, and CY190 (Callen et al., 1992i were screened to confirm the 
mapping of the STSs to the MEf-W interval within 16pl3,3. 

(ii) Microsatel lite- associated STSs: Cosmids in the minimal tiling 
path were screened for dinucleotide repeats to develop new genetic 
markers fronvthe region. In addition, 0163468 flizuka et al, 1993) 
and D16S261 7 were found to map to the MEFV region by recombi- 
nant analysis in FMF families (J. E. Balow ct al.. in preparation). 
DISS3275 had previously been isolated from iy23lE7 f26Fe7; French 
FMF Consortium. 1996) and was localized within cosmid 360H6 by 
oligonucleotide hybridization. 

Fluorescence in situ hybridization (FISH). Cosmids c383C12 and 
cRT70 were labeled with biotin- 14 -dATP (BioNick kit; Life Technol- 
ogies) and digoxigenin- 1 1-dUTP (Boehringer Manheim), respec- 
tively, using a nick-translation kit tLife Technologies). Two separate 
labeling reactions using both of the above-mentioned systems were 
performed on c57El, and the reactions were mixed prior to hybridiza- 
tion. EBV-Lransformed lymphoblastoid cell cultures from a normal 
male were treated with ethidivim bromide followed by colcemid and 
KGl treatment for mctaphase preparation. Dual-color FISH condi- 
tions were as described in Trask cf al. ) 1991) with minor modifica- 
tions; signal detection and amplification were done using the Signal 



Ainplificatinn Reagent Set (Oneor), Slides wiMe visualized '.i> 
Zeiss axiophot mictnscope with a dvial-pass filter. 

RESULTS 

Attempts to Construct a YAC Con tig of the FMF 
Candidate Region 

A number of CEPH megaYAC clones were previo 
mapped to 16pl3.3 near the vicinity MEFV ( Qhu 
kovetaL, 1995; Doggett e/' o/., 1995), Those were re 
rov D16S2844^ D16S2845, and 01683367 from cc 
(which contains the RFLP D16S246) and the c?: 
tetranucleotide repeat DI6S2622. In addition. 5 
were developed from cosmids that overlapped cR' 
on the centromeric side of MKFV (Table 1). Most o: 
YACs that mapped to 16pl3.3 were excluded from 
region because they were negative for these STSs f 
YACs, My655E10, My716Dl, My64lC2, and My64". 
were positive for D16S2622 and centromeric STSt 
negative for the c54G6 STSs. Only one YAC ci 
My806D3, was positive for markers on both ends oi 
interval. However, this YAC was found to be unsi 
since (i) 20 separate cultures each yielded differ 
sized YACs, ranging from 1 to 1.7 Mb and (ii) inter- 
PCR products of this YAC hybridized to far fewer t 
the expected number of chromosome-16 cosmids f 
YAC of expected size 1 to 1.7 Mb. The CEPH mega'-. 
library was also screened commercially with D 1 6S3 
and D16S3385 primers, which were derived i 
c54G6 and cRT4, respectively. No clones were idc 
fied with D16S3367 primers, while My64lC2 
My716Dl were selected by DI6S3385 primers Bot 
these YACs were later found to extend in the cer. 
nieric direction away from the MEFV locus. 

Due to the problems encountered with CE['H r::, 
YACs from (he region, we also screened smaller -n 
YAC libraries. The overall strategy we employed 
physical mapping is shown in Fig. 1. Briefly, cos- 
walking led to the development of new STSs. 'a ; 
were used to screen the ICI and Washington Univer 
YAC libraries. A total of 17 YAC clones were identi 
from these libraries (Fig. 2), with insert sizes rani 
from 150 to. 440 kb. Even though we did not corr.p 
a YAC contig of the region, these smaller YACs 
very useful as probes for the cosmid library scret?r. 
and in this manner contributed to the assembly of 
cosmid contig. 

Identification of BACs, PACs, and PI Clones 

A gridded BAC library became available only to^^ 
the end of the physical mapping project. BAC nl 
were screened with a 1.9-kb EcoRl fragment of c44 i 
and an SP6 (centromeric) end-clone probe from F 
273L24 to facilitate closure of the final gap in the 
mid contig (see below). Five positive clones were ide 
fied with the c441H9 fragment, one of which. F 
36P15, extended in the telomeric direction and 
eluded a region in which cosmids were found to cci^i 
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Starting Cosmids 



T3 and T7 Riboprobes 

I 

Gridded cosmid libraries 

I 

Identification of true positive clones 
by hybridization with fragments 
from starting cosmid 



Develop EcoRI restriction maps 
of positive clones 



Contig Assembly 



STS Primers from Genetic Markers 



Inter-alu PGR 
of positive clones 



Screening of 
YAC, PAC, BAG 
and PI libraries 



Riboprobe from 
the contig end 



Sequence of the 
contig end 



FIG. I. Strategy used in cloning and ptiysical mapping of the FMP" candidate region. The start.iiig cosmids were c54GR (contauiing r lie 
marker D16S246) and cRT70 (containing the markor D 1 682622 \- 



internal deletions. Using DNA from this BAG, we were 
able to estimate the size of this deleted region and to 
deduce its EcoKl restriction pattern. BAG 153 H7 was 
positive for the SP6 end of PAC 273L24 and extended 
just up to the deletion area ( Fig. 2), 

Ten PAC clones wore identified by screening with 
eight sets of primer pairs ( Fig. 2). PGR- and hybridiza- 
tion-based PI hbrary screening identified 15 clones 
(Fig. 2). Of these, Pi' B53E2 linked c54G6 to cosmids 
from contig 387 (which was known to contain the distal 
marker D16S94) and thereby helped orient the cosmid 
walk from this end. 

Construction of the Cosmid Contig 

A cosmid contig for the FMF region was generated 
by stepwise cosmid walking using riboprobes from the 
ends of cosmids and inter~A/u PGR probes from PI, 
PAC, BAG, and YAG clones, against both the finger- 
printed (1.5X) and the 16-2 (5.5x) chromosome-i6 cos- 
mid libraries. The initial attempts at cosmid walking 
were made in the fingerprinted library. This library 
provided 58% coverage of chromosome 16 in cosmid 
contigs and an additional 28% in fingerprinted single- 
ton clones (Staliings et al., 1992). The starting points 
for walking in this library were the cosmids containing 
the flanking genetic markers, i.e., c54G6 {D16S246) 
and cRT70 (D16S2622). At every step, overlaps were 
detected by hybridization of fragments from the start- 
ing cosmid and confirmed by detailed £^coRI and Notl 
restriction maps of each cosmid (Fig. 3). 

The first round of screening with cosmids c54G6 and 
cRTTO from each end of the interval was performed 
bidirectionally, using both T3 and T7 riboprobes, since 
the orientation of these clones was not known. We were 
able to orient the walk at cRT70 by finding overlap of 



T7-cRT70-positive clones with cRT8, which was 
known to extend in the centromeric direction (Petrij et 
al., 1995), and therefore walking was continued only 
from the T3 end of cRT70 toward the telomere. We 
were unable to identify any cosmids overlapping with 
c54G6 from the l.5x fingerprinted cosmid library. In 
subsequent walking experiments, gridded membranes 
of the 5.5 X library were also used. We were able to 
identify several overlapping cosmids with c54G6 from 
this library. The orientation of the walk from this end 
was established by showing that T3~c54G6-positive 
clones overlapped with PI B53E2, which had been de- 
termined to extend distally from c54G6. Thus, cosmids 
identified from the T7 end of c54G6 extended centro- 
meric across the FMF interval. 

As cosmid contig assembly was progressing from 
both ends of the candidate interva], D 1 6S468 was local- 
ized to the region by analysis of DNA from a somatic 
cell hybrid breakpoint mapping panel (Callen et al., 
1992, 1995). Analysis of recombinants placed this 
marker between D 16S246 and D16S2622, telomeric to 
MEF\^(J. E. Balowe/!oi., in preparation). Two cosmids, 
c363D9 and c420A3, were identified by screening cos- 
mid libraries with oligonucleotides for D16S468, and 
cosmid C28 {from which the D16S468 sequence was 
originally derived) was obtained from Dr. K. Hayashi. 
These three overlapping cosmids formed a small contig 
between the flanking markers for MEFV and were used 
as another nucleation point for cosmid walking. After 
extending this contig to c383G12 in one direction, and 
to c57El in the other direction, dual-color FISH analy- 
sis (Fig. 4) was performed to determine the relative 
order of these cosmids. Additional walks were initiated 
at cRT194 and cRT21 1 (cosmids provided by Dr. Breun- 
ing that were mapped by FISH telomeric of cRT70)- 
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In all, both chromosome-16 cosmid libraries were 
screened with 26 riboprobes; 14 inter-AZu PGR prod- 
ucts of PI, PAC, BAG, and YAC clones; one pair of 
oligonucleotides for the marker DI6S4G8; and DNA 
probes derived from the SP6 end of PAC 273L24 and 
the right (centromeric) end of iy23lE7 (26Fe7). These 
hybridizations identified 474 cosmid clones, of which 
179 mapped to the FMF candidate region, providing 
a continuous coverage of approximately 1050 kb. Only 
42 clones were analyzed from the fingerprinted li- 
brary, of which 23 mapped in the region. Twelve 
clones belonged to four previously characterized con- 
tigs, C211, C330, C664, and C387 (Doggett et at., 
1995). The remaining 11 cosmids identified in this 
library were singleton clones, nonoverlapping with 
cosmids in the previously available chromosome-16 
physical map. 

Based on the nature of the libraries used, we ex- 
pected to achieve six- to sevenfold coverage of the re- 
gion by cosmids. Although that turned out to be the 



case for most of the contig, areas of both underrepresen- 
tation and overrcpresentation were detected through- 
out the contig. It was particularly important to confirm 
overlaps in the underrepresented regions where small 
stretches of DNA at the ends of cosmids were minimally 
overlapping between only two cosmid clones (Fig. 3). 
These overlaps were confirmed by hybridization of end 
fragments to each other and to common fragments of 
overlapping Pi, PAC, and YAC clones. To estimate the 
size of the corresponding EcoRl fragment in genomic 
DNA, the degree of overlap between the end fragments 
was used for the overlap of the T3 end of c400Dl with 
the T3 end of c373C8. Restriction digests from YAG 
WyB34E5 were sized to determine the EcoRl fragment 
at the overlap of the T3 end of c363E3 with the T3 end 
of c363D9. Restriction digests from PAC clone 273L24 
were used to determine the EcoRl fragments at the 
overlap of the T7 end of c360A8 with the T7 end of 
c442E2/c442F3 and of the T7 end of c374H9 with the 
T3 end of c360H6. 
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Deletions and Instability in Cos m ids 

Although cosmids are relatively stable cloning vehi- 
cles, by detailed analysis of 179 clones encompassing 
this 1.05-Mb interval we found five cosmids from three 
separate regions with deletions. In one of these regions, 
two overlapping cosmids, c328C8 and c377A12, were 
unstable in culture. The deleted area in these cosmids 
is shown by a dashed line in Fig. 3. We were unable to 
identify any stable cosmid that spanned this region. By 
analyzing BAG 36P15 DNA, we estimated the extent 
of the deletion to be —30 kb, with EcoKl fragments of 
— 16, 7.5, 3.6, and 2.8 kb. Moreover, as is predicted 
from the location of Notl sites in cosmids that flank 
the deleted region, pulsed-field genomic Southern blots 
demonstrated an --llO-kh Notl fragment when hybrid- 
ized to a probe derived from c441H9 (data not shown), 
providing further evidence that BAG 36P15 contained 
an intact set of EcoRl fragments from this region. 

In a second region about 125 kb telomeric to this 
first unstable area, we found two overlapping cosmids 
(c360A8 and c373H9) that showed two different dele- 
tions relative to other cosmids in the region. In contrast 
to the c377A12/c328C8 segment, there are several sta- 
ble cosmids spanning this second region. The extent of 
the deletions in c360A8 and c373H9 was determined 
by hybridization of individual £^coRI fragments from 
c347C2 and c414G4. We estimated that -24 kb of geno- 
mic DNA had been deleted from c373H9 and 6.3 kb 
from C360A8. Multiple cultures of c360A8 and c373H9 
yielded reproducible restriction maps, indicating that 
the deletions observed in these cosmids occurred early 
in the library construction. 

Finally, in a region about 150 kb further telomeric, 
we found that cosmid G28 bears a deletion relative to 
other cosmids in the region. However, the significance 
of this observation is unclear, given the fact that this 
cosmid is derived from a different library, with a differ- 
ent donor of genomic DNA. 

Integration, of Physical and Genetic Maps 

In the course of constructing this map, we attempted 
to determine the physical locations of existing genetic 
landmarks and to identify new polymorphic markers 
from the region. At the telomeric end of the interval, 
the (AC),, repeat D16S94 fVK5) has been mapped to 
the \2-kbEcoRl fragment of c310G9, while the South- 
ern blot marker D16S246 (p218EP6) has been localized 
within c54G6. As is noted above, D26S468 maps to 
c363D9, and this information was used to establish a 



nucleation point for cosmid walking. We subsequei 
found that D16S3070 maps to the same EcoRl fr 
ment as does OI6S468, and sequence comparisons 
tween DWS468 primers and the D16S3070 amplii 
confirm that these are the same locus. Approxima^ 
150 kb centromeric to D16S468, CHLC microsatei 
D16S2617 was mapped to the 3.5-kb EcoKl fragrn 
of c414G4, D16S3275, which was recently reportec 
be the centromeric end of the candidate interval den 
by founder haplotypes (French FMF Consortiv 
1996), is approximately 75 kb further centrome 
within the 13.3-kb £:coRI fragment of c360H6. 

In the course of screening our clones for microsa 
Htes, we identified three polymorphisms that wore 
dependently identified by others. At the left (telomc 
end of iy58H4 f7Dh4), we identified a complex rep 
that, by sequence and map Jocation, is the same 
D16S475, In c49B4, about 20 kb centromeric 
DI6S468/DI6S3070, we identified an (AC), repea^ t 
is the same as Genethon marker D16S3082. Fin all 
tetranucleotide repeat polymorphism that we had 
veloped from cRT70 is the same as the CHLC mar 
D16S2622, 

In addition, we identified five microsatellites t 
have not been previously reported. These 
D16S3370, D16S3373, D16S3376, D16S3382, i 
D16S33S8. The first four have all been mapped wit 
our cosmid con tig, while DI6S3388 is derived from c 
approximately 700 kb centromeric to cRT70 fPetri 
aL, 1995). 

The boundaries of the FMF candidate interval I 
fined by non-Ashkenazi Jewish founder haplotypes 
DI6S4G8/D16S3070 and D 1 6S3275 (X^rcnch FMF C 
sortium, 1996). Historical recombinants in our pa 
of families slightly nariow the candidate interval to 
region between 01683082 and DI6S3373 (tf. E-. Ba 
et a/., in preparation). Depending on their exact p- 
tions within their respective EcoRl fragments, th 
markers are 190-205 kb apart. 

DISCUSSION 

In this article we present a clone-based, high-rest 
tion physical map of approximately 1 Mb of distal cl 
mosome 16p. The cloned interval includes the car 
date region for MEFV, the familial Mediterranean 
ver locus, and physically links extensive maps t 
have already been constructed distally in the polycy 
kidney disease (P/CD / )/tuberous sclerosis (TSC2) 
gion (Dackowski et aL, 1996) and proximally in 



FIG. 3. Detailed EcoBX and Notl restriction map of the ~1050-kb cosmid contig spanning the interval hetween D16S94 and DI6S2 
Cosmids are represented by horizontal lines with T3 and T7 ends marked on either side. In areas of overrepresentation, tiot all cost 
were restriction-mapped, and such cosmids are listed in boxes at tlie appropriate position in the contig. Genomic DNA is represented 
continuous horizontal line with its orientation relative to the chromosome indicated at the two ends. EcoRl sites are indicated by ver* 
lines going through the horizontal lines, whereas Notl sites are shown as lower vertical lines and an N. Numbers above the DNA hno 
EcoRl fragment sizes; numbers below are Notl-EcoHl or Notl-Notl fragment sizes. A shill (/) is used to list multiple EroRl siies if t 
relative order could not be determined. Deletions in cosmids are shown by dashed lines. Fragment sizes estimated from PAC27.TL24 
denoted as "6140". 
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Rubinstein-Taybi (/?57'S)/CREB-binding protein (CBP) 
region (Giles et ai, 1996). In toto, this completes ap- 
proximately 4 Mb of continuous sequence-ready maps 
of band 16pl3.3 CDoggett a/.. 1996). In the context of 
the chromosome-16 mouse/human somatic cell hybrid 
breakpoint map (Callen et ai, 1992, 1995), our 1.05- 
Mb contig completely spans the CY190-CY186/23HA 
mterval, with DJ6S94 (VK5) distal to the CY190- 
CY186 breakpoint and DI6S3381 and DI6S3383 proxi- 
mal to 23HA (data not shown). Presumably, normal 
sequences from both breakpoints are present in our 
contig. 

Within the 1.05-Mb cloned interval we have mapped 
27 STSs, including 18 markers not previously reported. 
This gives an average marker density of approximately 
1 STS every 40 kb, substantially better than the ge- 
nome-wide goal of 1 STS/100 kb. Moreover, within the 
region between D 16S468/D 1683070 and D16S475, the 
density is 1 STS every 23 kb. We have also developed 
an additional 8 STSs from the "^700-kb interval centro- 
meric to cRT70 (DI6S2622). 

Of the 27 STSs localized to the F'MF region cosmid 
contig, 10 are polymorphic microsatellites, and 4 
{DI6S3370, D16SS373, D16S3376, SindD16S3382] have 
not been described prior to the present report. Two of 
these latter microsatellites fall within the 250-kb FMF 
critical region recently defined on the basis of founder 
haplotypes {French FMF Consortium, 1996). We have 
also identified a fifth microsatellite (DI6S3388) from 
the more centromeric cosmid cN2. 

In addition, we found a second microsatellite from 
cN2 that, based on the sequence of the respective am- 
plimers, is the same locus as the Cienethon marker 
D16S3072 fAFMb015wa9). Thus, our directed search 
for microsatellites in this region identified a total of 
four polymorphisms {01633082, D16S475, D16S2622, 
and D16S3072) that had also been found through ge- 
nome-wide microsatellite mapping projects. Moreover, 
based on map location and primer sequence, we found 
that DI6S3070 (AFMa353yhl) is the same <x^Dl6S468 
(C28). Map location was important in raising our index 
of suspicion in this latter case, since amplimer se- 
quence is not available in the databases for D16S468. 
Our experience suggests that ascertainment of di-, tri-, 
and tetra nucleotide polymorphisms from this region 
may be relatively complete and emphasizes the utility 
of a detailed physical map in evaluating the uniqueness 
of newly identified markers. 

Based on CEPH family genotyping, the sex-averaged 
genetic distance between D16S246 and D16S468/ 
S3070 is ^1 cM, while the genetic distance between 
the latter locus and D16S475 is -2 cM (French FMF 
Consortium, 1996). On our map, the physical distance 
between D16S246 and D26S475 is approximately 635 
kb, indicating an approximate physical/genetic dis- 
tance ratio of 200 kb/cM. The whole chromosome aver- 
age for chromosome 16 is 95 Mb/152 cM or 600 kb/cM 
(Kozman et aL, 1995). It therefore appears that the 
FMF region is a "hot spot" of recombination relative to 



the rest of chromosome 16 and markedly so relative 
to the genome-wide average of 1000 kb/cM. Since the 
physical distance from D16S246 to O16S468/S3070 is 
roughly half the physical distance from 0I6S468/ 
S3 070 to D16S475, the physicaUgenetic distance ratio 
is relatively uniform over the interval, to the present 
level of resolution. 

In addition to being relatively recombinogenic in 
man, this 1-Mb region is relatively unstable as large- 
insert YACs. Similarly, Doggett et aL (1995) have re- 
ported an apparent lack of megaY'AC clones in the dis- 
tal 2.85 Mb of 16pl3.3 (telomeric to MEFV), while in- 
stability of YACs in proximal 16pl3.3 (centromeric to 
iWEP^l led Petrij et ai (1995) to use cos m ids in their 
hunt for the Rubinstein-Taybi gene. Taken together, 
these data imply that the distal --4 Mb of 16p is rela- 
tively unstabie/underrepresented in megaYACs. We 
were, however, able to identify stable smaller-insert 
YACs from the I CI library in both the proximal and 
the distal regions of our contig, while the central region 
of the contig, between iy231E7 (26Fe7) and iy237E7 
f27Ce7), was not represented. It is intriguing to note 
that we also observed instability and relative underrep- 
resentation of cosmids in this same central region. 
Thus, there may be sequences in this particular region 
that i mpair stable propagation in both yeast and bacte- 
rial systems. The failure to find any stable megaYACs 
spanning the interval may be due to the fact that any 
such clone would necessarily contain these sequences. 
Based on data from tlie PKDl region, one can speculate 
that instability may be due to the presence of repeats 
and duplications and/or a high CC content. It is un- 
likely that the instability in YACs is directly related to 
the high recombination frequency observed across this 
interval, since another hot spot at l6pL2 is well repre- 
sented in megaYACs (Callen et a!., 1996). 

The utility of BAC clones in spanning this unstable 
region suggests that low copy number may be im- 
portant in propagating these sequences in bacterial 
systems. The close agreement of physical distance esti- 
mates based on BAC clones and pulsed-field mapping 
of genomic DNA argues against any major additional 
gaps in our map, In addition, this experience once again 
underscores the importance of employing several com- 
plementary cloning systems in constructing an accu- 
rate map of this magnitude. 

Of course, the rat son d'etre for this map is the identi- 
fication, by positional cloning, of the gene causing FMF. 
We are currently employing exon-trapping, direct 
cDNA selection, and single-pass sequencing to develop 
a detailed transcriptional map of the region. Given that 
distal 16p is relatively gene dense, it is likely that the 
resources created here will be of use in other cloning 
projects, as well. 
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i\'o(i' added tn proof- Upon analysis of the sc<]ucncc general od frnni 
(;'J99Ain, wo found an Fa:o\\\ fraf^ment of 3.1 kb that is i^ot shown 
in Fig. 3. This fragment lies between the 3.8- and the 3-kh fragments 
ai)oul 15 kb telomeric to the marker f)!6S3370. The fraKniont was 
not deteeted by ori^jinal arialysis <lue to the presence of two <il.hor 
fra^^nients of the same size in this I'effion. 
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